Stable monopole-antimonopole string background in SU{2) QCD 
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Motivated by the instability of the Savvidy-Nielsen-Olesen (SNO) vacuum we make a systematic 
search for a stable magnetic background in pure 5(7(2) QCD. It is shown that a pair of axially 
symmetric monopole and antimonopole strings is stable, provided that the distance between the 
two strings is less than a critical value. The existence of a stable monopole-antimonopole string 
background strongly supports that a magnetic condensation of monopole-antimonopole pairs can 
generate a dynamical symmetry breaking, and thus the magnetic confinement of color in QCD. 
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1. Introduction 

It has long been argued that the monopole conden- 
sation could explain the confinement of color through 
the dual Meissner effect Jj. Indeed, if one assumes the 
monopole condensation, one could easily argue that the 
ensuing dual Meissner effect guarantees the confinement 
0. There have been many attempts to prove this sce- 
nario in QCD Unfortunately, the earlier attempts 
had encountered the well known problem of Savvidy- 
Nielsen-Olesen (SNO) vacuum instability j^. In fact the 
effective action of QCD obtained with the SNO vacuum 
develops an imaginary part, which implies that the vac- 
uum is unstable 0, IH . This instability of the magnetic 
condensation has been widely accepted and never been 
convincingly revoked. 

Recently it has been shown that, if one uses a proper 
infra-red regularization which respects causality, the 
imaginary part in the effective action can be removed 
0- In addition, a numerical evidence for the stable mag- 
netic condensation has been found in lattice simulation 



The purpose of this paper is to search for a sta- 
ble classical magnetic background in SU{2) model of 
QCD. We analyze the stability of two classical mag- 
netic backgrounds, a pair of axially symmetric monopole- 
antimonopole strings and a pair of magnetic vortex- 
antivortex strings, and show that the first one becomes 
stable provided the distance between two strings is less 
than a critical value. As far as we understand, the pair of 
axially symmetric monopole-antimonopole strings consti- 
tutes a first explicit example of a stable magnetic back- 
ground in QCD. More importantly, the result can serve as 
a strong indication that a monopole-antimonopole con- 
densation can provide a stable vacuum in QCD. 
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2. Instability of Wu-Yang monopole back- 
ground 

Let us start with a gauge invariant Abelian projection 
in SU{2) model of QCD which includes explicitly a topo- 
logical degree of freedom expressed by the unit isotriplet 
n. We decompose the gauge potential into the restricted 
potential and the off-diagonal part (valence gluon) 

X, 



1 

9 

1, fi ■ X,, 



—h X df^h 



Xfj. — A^ + Xfj,, 



0), 



(1) 



Notice that the restricted potential is precisely the con- 
nection which leaves fi invariant under the parallel trans- 
port. 



D^n = df^n -(- gA^ x n = 0. 



(2) 



The restricted potential A^ has a dual structure which 
can be seen from the field strength decomposition 



F^u = dfj^Au — dpAp^, 

where C*^ is the "magnetic" potential 0- 
With the decomposition one has 

Ffj.1, = F^u + Dfj^Xi, — DijX^ + gXf^ x X^, 
so that the Lagrangian can be written as follows 
1 



(3) 



(4) 



c = --^Fl,--^{b^x,-bat.f 



2 ^" 



{X^ X X,) -^-{X^x X,f 



(5) 
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With the gauge fixing condition D^X^ = the one-loop 
correction to the effective action reduces to the form 

exp(iA5) = Viei^'^K^,^ Det^Mpp, 

Mfp = DD, = + ig{Af, + (7^), (6) 

where the operators K^i, and Mpp originate from the 
functional integration over the off-diagonal gluon and 
Faddeev- Popov ghost respectively (the contribution from 
integration over the quantum part of the Abelian field 
Af^ is trivial). For arbitrary static magnetic background 
Ffj^i, + H^u one can simplify the one-loop correction 

AS" = ilnDet(-i)2 + 2a) + i lnDet(-I)2 - 2a), (7) 



where a = g\j-^H^^, hereafter, for the brevity of nota- 



tion, we employ a single notation H^i, for the additive 
combination F^^, + H^^. 

Before we discuss the stability of monopole-antimon- 
opole pair, we first review the instability of the Wu-Yanj 
monopole because two problems are closely related 
The Wu-Y ang mo nopole solution of charge q/g is de- 
scribed by 



ang mc 



An = fi X 

9 

(sin 9 cos q(j) \ 
sin 9 sin q(f) 
cos 9 J 



(8) 



here, a = is given by the magnetic field strength of the 

Wu-Yang monopole. The absence or presence of negative 
modes of the operator K implies stability or instability 
of the classical background against small fluctuations of 
the gauge potential. To calculate the eigenvalues of the 
operator K one can write down the eigenvalue equation 
as the following Schroedinger type equation for a complex 
scalar field which plays a role of wave function 



K^{r,9,(j)) = E-^{r,9,(j}), 
iq 



K = -A -2 



1 



2 

, . 2n COS Od^ + % cot^ 61 ± 2-1, 
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L2 



-deism 9 de) + 
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sin^ I 



-^dr{r dr) ^, 



(11) 



where L is the angular momentum operator. Notice that 
here the ± signatures represent two spin orientations of 
the valence gluon. 



Substituting the separated solution 
■^[r,9,<P)^R{r)Y{9,<i>), 



(12) 



into 111|) one obtains the equation for the angular eigen- 
function Y{9, (f>) 



2iqcoa9 



sm 



d^ + q^ cotH)Y{9, 



^\Y{9,4>). 



(13) 



Moreover, with 



where (r, 9, 4>) is the spherical coordinates and q is an 
integer monopole charge. In the Abelian formalism it 
is more convenient to describe the monopole in terms 
of the magnetic potential C^, This implies the compo- 
nents of the magnetic potential and the magnetic field 
strength H to be as follows 



-t-oc 



Cu. 



H, 



(cos 6* - l)d^(j), 



gr 

q x"' 



(9) 



We use the parametrization for the magnetic potential 
slightly different from the parametrization in 
where the magnetic potential is defined on two coordi- 
nate patches. 

To study the stability of the monopole background 
we should consider an operator obtained by taking the 
second variation of the classical Lagrangian with respect 
to small fluctuations of the field X^. The operator is 
identical to the operator i^^jy, ©, and the problem of 
flnding unstable modes is reduced to calculation of the 
scalar functional determinants in 



DetT^ = Det(-i)2 ± 2a), 



(10) 



*«(<?!>) = - 
one can reduce (|13|l to 

^ ^ -de{sm9de) 



/2n 



exp(im0). 



(14) 



sm( 



(m + q cos ( 
sm^9 



e = Ae. (15) 



This is exactly the eigenvalue equation for the monopole 
harmonics which has been well-studied in the literature 
[Tol IT^ . From the equation one obtains the following 
expression for the monopole harmonics and the corre- 
sponding eigenvalue spectrum 



r,,„((?,0) = e,,„,(^)$,^(0), 

|m ± q 



qjm{0) = il-COs9)''+{l+COs9)''-Pk{cOS( 



7± 

A=j(i + l)-9' 

j = k + 'y+ +7_, 



fc = 0,l,2,... 



(16) 



where Pk{x) is the Legendre polynomial of order k. The 
quantum number j is analogous to the orbital angular 
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momentum quantum number / of the standard spherical 
harmonics Yi^ , except that here j starts from a non-zero 
integer value for a non- vanishing monopole charge q. 

Now, consider the equation for the radial eigenfunc- 



tion 



^T^r'^) + 1) - q'] T ^+E)Rir) 

ar ar r / 



0. 



(17) 



With R{r) — -x(r) one obtains 



{-^-^ + 1) - 9' ± 2?] + E)x{r) - 0. (18) 

The equation has a general solution in terms of Bessel 
functions of the first kind Ju{z) 



V = -^l + A[j{j + l)-q^±2q], (19) 

where Ci [i — 1, 2) are the integration constants. For 
positive values of v and E the finite solutions oscillat- 
ing at the infinity and vanishing at the origin are given 
by Ci = 0. The negative eigenvalues of E can come 
only from (|18|l with the lower negative sign (which cor- 
responds to the operator —D"^ — 2g/r^) and the lowest 
value of J = 1 with q = 1. In this case we have to solve 
the equation 



E x = 0, 



(20) 



which is nothing but the one-dimensional Schroedinger 
equation with the attractive potential — 1/r^ ^M- The 
solution to this equation has a continuous eigenvalue 
spectrum for both positive and negative energies, and 
leads to the radial eigenfunction R{r) which behaves like 



const 




(21) 



near the origin. The solution has an infinite number of 
zeros approaching the point r = 0. This can be inter- 
preted as a valence gluon movin g ar ound the monopole 
and falling down to the center |3, UM ■ 

One can observe that for the lowest energy states, 
(j — 1), the undesired attractive potential proportional 
to — 1/r^ in H18|l vanishes when q = 0. This can serve 
as a hint that one might expect the absence of negative 
modes for a magnetic background with zero monopole 
charge. The simplest magnetic configuration with a to- 
tal vanishing monopole charge can be realized as a Wu- 
Yang monopole-antimonopole pair. Unfortunately, since 
near the location of the (anti-)monopole we still have the 
attractive potential part — 1/r^ one can verify that the 



Wu-Yang monopole-antimonopole pair has to be unsta- 
ble. 

3. Axially symmetric monopole string 

The axially symmetric monopole string can be re- 
garded as an infinite string carrying homogeneous 
monopole charge density along the string. The magnetic 
field strength of the axially symmetric monopole string 
can be written in a simple form in the cylindrical coor- 
dinates (/9, 0, z) 



C,, 



0, 



A,. 



H = -p, 

p 



(22) 



where a is the monopole charge density and r is an arbi- 
trary constant which represents the translational invari- 
ance of the magnetic field along the z-axis. Just like the 
monopole solution ijHJl the above monopole string forms a 
singular classical solution of SU{2) QCD. In the follow- 
ing we will assume a = 1 and r = for simplicity, since 
this will not affect the stability analysis. 

Let us consider the eigenvalue problem for the opera- 
tor K 



K^{p,4),z)^E^{p,4>,z). 



With 



+ 00 



H Fm{p,z)<^ra{(t>): 



(23) 
(24) 



ra— — oo 



and repeating the steps of the previous section we obtain 
the following eigenvalue equation 



pp 



En ^" F^^ 



[m — z)^ 



± - -E 
P 



F = Q. (25) 



By shifting z to z + m one can put m = 0. The quantum 
mechanical potential of this equation behaves like ±2/p 
near p = 0. So we still have an undesired attractive po- 
tential —2/ p. This implies two things. First, the attrac- 
tive interaction of the axially symmetric monopole string 
background is less severe than the attractive interaction 
of the spherically symmetric monopole background. So 
we can expect the absence of continuous negative en- 
ergy spectrum for the axially symmetric monopole string 
background. Secondly, the attractive potential —2/p tells 
that the monopole string background must still be unsta- 
ble, because it indicates the existence of discrete bound 
states with negative energy. 

To confirm this we make a qualitative estimate of the 
negative energy eigenvalues of 1)25(1 . We look for a solu- 
tion which has the form 



E{P,Z) = ^ fn{p)Zn{a 



VP' 



(26) 
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where Hn{x) is the Hermite polynomial. Notice that 
Z„{x) forms a complete set of eigenfunctions of the har- 
monic oscillator, 

(^-a;2)^„(a:) = -(2n+l)Z„(a;). (27) 
Substituting ^ into ^ we obtain 

oo 

{( 



Cr/n ^ /I ^ f2)^W ^dfndKn 

^-^ ^ V dp^ ^ p p^' dp / " dp dx 



■^"Up3 dx^ ^p'^Vp P dx J 



2n + 1 ± 2 

~3 



V4p4 4p3^ p" p 
= 



E 



(28) 



Using the recurrence relations and orthogonality proper- 
ties of Hermite polynomials one can derive the equations 
for fnip) 



f d^ 2n + Z d 4^2 - 2n - 1 

\ ^ 

V dp'^ 2p dp 

2n + 1 ± 2 



16p2 



P 



E f, 



r/n — 4 + „o 9 fn- 



l / d n-l 



64p2''" " 32^2-"'-" 4pVdp 4p 
3_ , 3(n + l)% 



/n- 



16p 

(n-M)(7i + 2) / d 3n + 2\ 
(n+ l)(n + 2)(n-f 3) 



4p2 



where /« = for negative integer n. So we have infi- 
nite number of equations for infinite number of unknown 
functions fn{p)- 

Notice that the left hand side of the last equation con- 
tains a second order differential operator with the quan- 
tum mechanical potential 



U ^ 



2n + 1 ± 2 



(30) 



The potential becomes attractive only if n = 0, so that 
in a first approximation we expect that the negative en- 
ergy eigenvalues will originate mainly from the lowest 
bound state with n = of the harmonic oscillator part. 
So, by neglecting all /„ with n ^ we can still get an 
approximate qualitative solution for /q. In such an ap- 
proximation the equation reduces to a simple one 



+ 



3 d 1 



1 



dp^ 2p dp p 16p2 



+ S/o = 0. (31) 



The solution to this equation has a new integer quantum 
number k, 



foAp) 



Eh — 



ai+i 




2s + 3/2)2' 

y/\E^\{2l + 2s + 3/2) - 
(/ + l)(/ + 2s + 3/2) 



-ai. 



(32) 



With this we may express the corresponding eigenfunc- 
tion ^'fc as 

z2 

*fc(p,0,z) = iVfcexp( - — )/o,fc(p), (33) 

where Nk is a normalization constant. One can find the 
lowest energy eigenvalues 



Eo 
El 
E2 
E3 
Ei 



-0.343. 
-0.073. 
-0.031. 
-0.017. 
-0.011. 



(34) 



This confirms that the axially symmetric monopole string 
background is indeed unstable. 

Surprisingly, we find that the approximate solution 
(|33|1 can also be obtained as an exact solution of varia- 
tional method with the trial function F of the form 



(/n+3-3(n + 4)/„+4), (29) p^p^^^ 



= Np' exp (- f3kP' 



l = k 

^aip\ 

1=0 



(35) 



where s,/3fc,7,a; are treated as variational parameters. 
In other words, the variational minimum of the energy 
functional with the above trial function is provided ex- 
actly by the solution 

To quantify the accuracy of our approximate analytic 
solution we solve numerically the starting equation H25(l 
(with the lower negative sign and m = 0). The obtained 
numerical solution for ^E* have the same essential singular- 
ity structure as in H33|l . The corresponding lowest energy 
eigenvalues 



Eo = 


-0.545 


El = 


-0.093 


E2 = 


-0.036 


E3 = 


-0.019 


Ei^ 


-0.011 



(36) 

confirm that the solution H33|) provides a good qualitative 
estimation of the energy spectrum to analyse the vacuum 
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stability of the axially symmetric monopole string. What 
is more important is that we can apply the structure of 
that solution to the analysis of the stability problem for 
the monopole-antimonopolc string configuration in the 
subsequent section. 

4. A stable magnetic background 

The main idea how to construct a stable magnetic 
background is quite clear. Consider a pair of axially sym- 
metric monopole and anti-monopole strings which are or- 
thogonal to the xy-plane and separated by a distance a. 
Due to the opposite directions of the magnetic fields of 
the monopole and anti-monopole strings the attractive 
part of the quantum mechanical potential U{p) in the 
eigenvalue equation falls down as U{p) 0{—a/p^) 
when p — > cxD. This allows the centrifugal potential to 
be competitive with the attractive part and prevail for 
small enough values of a. By decreasing the distance a 
we can decrease the effective size of the quantum me- 
chanical potential well, so that the bound state energy 
levels will be pushed out from the well. This, with the 
positive asymptotics of the potential at infinity, implies 
that the bound states will have disappeared completely 
at some finite critical value of a. 

To show this, consider a pair of axially symmetric 
monopole and anti-monopole strings located at {p — 
a/2,0 = 0) and (p = a/2,0 = tt) in cylindrical coordi- 
nates. The magnetic field strengths H± for the monopole 
and anti-monopole strings are defined as follows 



form 



P± = 



a 

± — P±, 
P± 

a 

P± 
^ 2' 



p ± ap cos ( 



4 ' 



(37) 



where a is the two-dimensional vector starting from the 
anti-monopole string to the monopole string in xy-plane. 
From now on we will assume a — 1 without loss of gen- 
erality. The total magnetic field is given by 



H = H. 



Hp 
H 



p cos < 

' 2 



a 

p + — cos ( 



P 



P+P- 



(38) 



One can express the corresponding vector potential in 
terms of -ffp,0 components 

Af, = zH^d^p - pzHpd^,4>. (39) 

The eigenvalue equation for the operator K takes the 



{ - -dp{pdp)--dl~dl~2i'-(H^dp~Hpd^^ 

z^H^ ±2H^^{p,(l},z) = E^{p,(l3,z). (40) 



The equation can be interpreted as a Schroedinger 
equation for a massless gluon in the magnetic field of 
monopole and anti-monopole string pair. 

Let us analyse the equation qualitatively first. We 
will concentrate on the potentially dangerous term —2H 
in (|40|l . The singularities of the term z^H^ determine the 
essential singularities of the differential equation. One 
can extract the leading factor of the solution and look for 
a finite solution for the ground state in the form similar 
to 



■^{p,(j>,z) = (np+p^) 3exp(---^ )F{p,(j)), (41) 

^P+P- 



where F{p, (jj) is normalized by 



\Fip, ^)\^pdpd^ = l. 



(42) 



The solution describes a wave function localized mainly 
near the string pair. The wave function vanishes exactly 
on the axes of the strings. This implies that the ground 
state has a non-zero orbital angular momentum which 
provides a centrifugal potential as we will see later. 

The lowest negative eigenvalue of this equation can be 
obtained by variational method by minimizing the corre- 
sponding energy functional 



E 



1 



1 



dp{pdp)--di-di-2z-iH^dp 



-Hpd^) + z^H^ ± 2H ^dpdzd4>. 



Now, with 



-t-oo 



F(p,0)= ^/,„(p)$„(0), 



(43) 



(44) 



one may suppose that the main contribution to the 
ground state energy comes from the first term of Fourier 
expansion with to = 0. With this we can perform the 
integration over z-coordinate and simplify the above ex- 
pression to 



Pfip) {-dpp- -dp + C/(p, 0)) f{p)dpdc^, 



U{pA)^ 



2ap+p- 



2pipi 



(45) 



where /(p) = /o(p) and C/(p, 0) is an effective potential. 
Since the energy eigenvalues decrease with decreasing the 
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parameter a, to study the features of the potential at 
small a we make the following rescaling 



P - 
E 



ap, 



f^f/o 



E/a" 



(46) 



Under this rescaling the potential near the origin can be 
approximated to 



U{p,(j)) -4a + (8- 16acos20)p2 



(47) 



So that the potential reduces to a two-dimensional har- 
monic oscillator potential whose depth decreases as a 
goes to zero. This implies that the negative energy eigen- 
values can disappear for some small values of a if the 
asymptotics of the potential at infinity becomes positive. 
To see this we perform the integration over the angle 
variable 4> in the energy functional, and with the change 
of variable 



(48) 



we obtain the following equation which minimizes the 
energy, 



dp2 



V{p) x{p) 

1 



Ex{p). 
8p2 



V(a4 - 16/)2 
8a 16aV^ 



7r^(a2 - 4/92): 



(a2 -4p2)2' 



(49) 



where K{x) is the complete elliptic integral of the first 
kind. Notice that the first term in the potential does not 

produce bound states because the potential leads to 

P 

negative energy eigenvalues only for n > 1/4. From the 
last equation we obtain the asymptotic behavior of the 
potential V{p) near space infinity 



V{p) 



4 



(50) 



This tells that the potential becomes positive when the 
distance a becomes less than the critical value acr (in the 
unit 1/a) 



a < Cc 



(51) 



To check the analytic estimate of the critical value acr 
we solve numerically the original equation (|in|) . Since the 
wave function ^{p, (j), z) becomes spread to long distance 
for small energy eigenvalues we extend sufficiently the 
upper limit of p in the domain (0 < p < Pupper,0 < 
(j) < 2tt) and increase the mesh near the location of the 
monopole and antimonopole strings. The extrapolation 
of the results from finite small energy eigenvalues to zero 
gives the following critical value within 6% of accuracy 



0.246. 



(52) 



The numerical result confirms that qualitatively the ap- 
proximate solution (|41|l describes the correct physical 
picture. In particular, this tells that a pair of monopole 
and antimonopole strings becomes a stable magnetic 
background if the distance between two strings is small 
enough. 

5. On stability of vortex-antivortex pair 

Recently an alternative mechanism of confinement has 
been proposed which advocates the condensation of mag- 
netic vortices However, it has been known that 
the magnetic vortex configuration is unstable IT5l |. So it 
would be interesting to study the stability of the vortex- 
antivortex pair. In this section we study the stability of 
a special vortex-antivortex configuration. 

Let us start with a single magnetic vortex given by 

1 , 



H 



A,. 



pd^,(, 



(53) 



Notice that, unlike the monopole string (|22(l . the vor- 
tex configuration is not a classical solution of the sys- 
tem. But since such a type of configuration multiplied 
by an appropriate profile function has been studied in the 
framework of various approaches, we will consider that 
special vortex configuration in the following. 

One can write down the corresponding eigenvalue 
equation of the operator K 



92 



1 



d„ 



1 



2i 



£;*(p,0,z). 



*(p,0,2) 

(54) 



The equation becomes separable in all three variables. 
With factorization 



*= 5]/™(p)5(^)'i>m(0), 9{z) = l, (55) 

one obtains the following ordinary differential equation 
for f{p) from 



m 2 . 2 



{i + -y±--E fip) = o. 
p p 



(56) 



The bound states are possible for the potential — 2/p with 
non-positive integer m, in which case the corresponding 
solution can be obtained 



/n,™(p) = pl"le- 



ApI 



Ap) 



En.m k 
flfe P > 

fc=0 



y/l - En,m{2k + 2|to| + 1) - 2 -|- 2m 



E,, 



0,1,2, 



(fc + l)(fc- 
4(1 -m)2 
(2n-|-2|m|-|-l)2' 
..; m = 0,-l,-2, 



2|m| + l) 



(57) 



7 



Clearly, the ground state has a negative energy Eq^q, 
which tells that the vortex configuration is unstable. 

There is a principal difference between the axially 
symmetric monopole string and the vortex configuration. 
The ground state of the monopole string has a non-trivial 
centrifugal potential. But the ground state eigenfunction 
/o,o(p) of the vortex configuration corresponds to an S- 
state, which implies the absence of the centrifugal poten- 
tial. This plays a crucial role in the existence of the neg- 
ative energy eigenstates in the case of vortex-antivortex 
pair. 

The vortex-antivortex background is described in a 
similar manner as the monopole-antimonopole string 
background in the previous section. The potential has 
the form 



H 



1 



1 



-ip- 



■ cos ( 



-ip- 



■ cos ( 



1 

p- 



(58) 



where a is the distance between the axes of the vortex 
and anti-vortex. The eigenvalue equation corresponding 
to the operator K is partially factorizable in z-coordinate 



+Al ± 2H 

h'' 



dl- 



2i{Apdp + -^A^d^) 



(59) 



The numerical analysis of the equation shows that 



there is no critical value for the parameter a, so that 
the negative energy eigenvalues exist for any small a. 
Qualitatively one can see this from the effective potential 



V - At 



2H. After averaging over the angle variable 



one can find the asymptotic expansion of the potential 
near the origin and infinity 



Stt 



V{p) 



64 

—P- 



IGtt 



^ 2 ' 



-^P 



{p ~ oo). 



(60) 



This shows that there is no centrifugal term which could 
prevent the appearance of bound states for small a. 
Whether the instability problem can be overcome with 
a more complicate configuration of the vortex-antivortex 
remains an open question. 

In conclusion, we have shown that the axially sym- 
metric monopole-antimonopole string background is sta- 
ble under the small field fiuctuation if the distance be- 
tween two strings becomes less than the critical value 
flcr — 1/4. The existence of the stable classical mag- 
netic background implies that "a spaghetti of monopole- 
antimonopole string pairs" could generate a stable vac- 
uum condensation. This would allow a magnetic confine- 
ment of color in QCD. 
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